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Abstract 

A measurement of the effective B® — > K + K~ lifetime is presented using approxi- 
mately 37 pb" 1 of data collected by LHCb during 2010. This quantity can be used 
to put constraints on contributions from processes beyond the Standard Model in 
the Bg meson system and is determined by two complementary approaches as 

tkk = 1-440 ± 0.096 (stat) ± 0.008 (syst) ± 0.003 (model) ps. 
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1 Introduction 



The study of charmless B meson decays of the form B — > h + h'~ , where is either a 
kaon, pion or proton, offers a rich opportunity to explore the phase structure of the CKM 
matrix and to search for manifestations of physics beyond the Standard Model. The 
effective lifetime, defined as the decay time expectation value, of the B° s meson measured 
in the decay channel B® — > K + K~ (charge conjugate modes are implied throughout the 
paper) is of considerable interest as it can be used to put constraints on contributions 
from new physical phenomena to the B® meson system [l}j4]. The B® — > K + K~ decay 
was first observed by CDF |5,6|. The decay has subsequently been confirmed by Belle [7]. 

The detailed formalism of the effective lifetime in B° s — > K + K~ decay can be found in 
Refs. (3j|4]. The untagged decay time distribution can be written as 

V{t) oc (l-^ A rJe- r ^ + (l + ^ Ar Je- r ^. (1) 

The parameter A&r s is defined as A A r a = — 2Re(A)/ (1 + |A| 2 ) where A = (q/p)(A/A) 
and the complex coefficients p and q define the mass eigenstates of the B®-B° s system 
in terms of the flavour eigenstates (see, e.g., Ref. pj), while A (A) gives the amplitude 
for B® (B° s ) decay to the CP even K + K~ final state. In the absence of CP violation, 
Re(A) = 1 and Im(A) = 0, so that the distribution involves only the term containing T^. 
Any deviation from a pure single exponential with decay constant T^ 1 is a measure of CP 
violation. 

When modelling the decay time distribution shown in Eq. [T] with a single exponential 
function in a maximum likelihood fit, it converges to the effective lifetime given in Eq.[2]|9|. 
For small values of the relative width difference Ar s /r s = (T L — T H )/ {{T L + r#)/2), the 
distribution can be approximated by Taylor expansion as shown in the second part of the 
equation [3] 



t~kk = t b o 



i-y. 
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1 + 2A Ars y s + yj 
1 + A A r s y s 



r B o(l+A Ars y s + 0(y 2 s )), (2) 



where r^o = 2/ (r# + Tl) = T s 1 and y s = Ar s /2r s . The Standard Model predictions for 
these parameters are „4 A r s = -0.97^;^ § and y s = 0.066 ± 0.016 

The decay B® — > K + K~ is dominated by loop diagrams carrying, in the Standard 
Model, the same phase as the B° s -B° s mixing amplitude and hence the measured effective 
lifetime is expected to be close to T^ 1 . The tree contribution to the B®— > K + K~ decay 
amplitude, however, introduces CP violation effects. The Standard Model prediction is 
t kk = 1-390 ± 0.032 ps 13). In the presence of physics beyond the Standard Model, 
deviations of the measured value from this prediction are possible. 

The measurement has been performed using a data sample corresponding to an inte- 
grated luminosity of 37 pb~ collected by LHCb at an energy of \fs = 7 TeV during 2010. 
A key aspect of the analysis is the correction of lifetime biasing effects, referred to as the 
acceptance, which are introduced by the selection criteria to enrich the B meson sample. 
Two complementary data-driven approaches have been developed to compensate for this 
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bias. One method relies on extracting the acceptance function from data, and then ap- 
plies this acceptance correction to obtain a measurement of the B° s — > K + K~ lifetime. 
The other approach cancels the acceptance bias by taking the ratio of the B° s — > K + K~ 
lifetime distribution with that of B°— > K + n~ . 

2 Data Sample 

The LHCb detector III] is a single arm spectrometer with a pseudorapidity acceptance 
of 2 < f] < 5 for charged particles. The detector includes a high precision tracking 
system which consists of a silicon vertex detector and several dedicated tracking planes 
with silicon microstrip detectors (Inner Tracker) covering the region with high charged 
particle multiplicity and straw tube detectors (Outer Tracker) for the region with lower 
occupancy. The Inner and Outer trackers are placed after the dipole magnet to allow the 
measurement of the charged particles' momenta as they traverse the detector. Excellent 
particle identification capabilities are provided by two ring imaging Cherenkov detectors 
which allow charged pions, kaons and protons to be distinguished from each other in the 
momentum range 2-100 GeV/c. The experiment employs a multi-level trigger to reduce 
the readout rate and enhance signal purity: a hardware trigger based on the measurement 
of the energy deposited in the calorimeter cells and the momentum transverse to the 
beamline of muon candidates, as well as a software trigger which allows the reconstruction 
of the full event information. 

B mesons are produced with an average momentum of around 100 GeV/c and have 
decay vertices displaced from the primary interaction vertex. Background particles tend 
to have low momentum and tend to originate from the primary pp collision. These features 
are exploited in the event selection. In the absolute lifetime measurement the final event 
selection is designed to be more stringent than the trigger requirements, as this simplifies 
the calculation of the candidate's acceptance function. The tracks associated with the 
final state particles of the B meson decay are required to have a good track fit quality 
(x 2 /ndf < 3 for one of the two tracks and x 2 / n df < 4 for the other), have high momentum 
(p > 13.5 GeV/c), and at least one particle must have a transverse momentum of more 
than 2.5 GeV/c. The primary proton-proton interaction vertex (or vertices in case of 
multiple interactions) of the event is fitted from the reconstructed charged particles. The 
reconstructed trajectory of at least one of the final state particles is required to have a 
distance of closest approach to all primary vertices of at least 0.25 mm. 

The B meson candidate is obtained by reconstructing the vertex formed by the two- 
particle final state. The B meson transverse momentum is required to be greater than 0.9 
GeV/ c and the distance of the decay vertex to the closest primary pp interaction vertex 
has to be larger than 2.4 mm. In the final stage of the selection the modes B®— > K + K~ 
and B° — > K + tt~ are separated by pion/kaon likelihood variables which use information 
obtained from the ring imaging Cherenkov detectors. 

The event selection used in the relative lifetime analysis is very similar. However, 
some selection criteria can be slightly relaxed as the analysis does not depend on the 
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Figure 1: Results of the relative lifetime fit. Left: Fit to the time- integrated KK mass 
spectrum. Right: Fit to the KK decay time distribution. The black points show the 
total number of candidates per picosecond in each decay time bin, the stacked histogram 
shows the B®—> K + K~ yield in red (dark) and the background yield in grey (light). 

exact trigger requirements. 



3 Relative Lifetime Measurement 

This analysis exploits the fact that the kinematic properties of the B° s — > K + K~ decay 
are very similar to those of B°— > K + n~ . The two different decay modes can be separated 
using information from the ring imaging Cherenkov detectors. The left part of Fig. [T] 
shows the invariant mass distribution of the B° s — > K + K~ candidates after the final 
event selection. In addition 1,424 B° — > K + n~ candidates are selected. Using a data- 
driven particle identification calibration method described in the systematics section, the 
remaining contamination in the B®— > K + K~ sample from other B— > h + h'~ final states 
in the analysed mass region is estimated to be 3.8%. 

B mesons in either channel can be selected using identical kinematic constraints and 
hence their decay time acceptance functions are almost identical. Therefore the effects of 
the decay time acceptance cancel in the ratio and the effective B° s ^t K + K~ lifetime can 
be extracted relative to the B°-> K + n~ mode from the variation of the ratio R(t) of the 
yield of B meson candidates in both decay modes with decay time : 

R(t) = J R(0)e" i ( T ^ - T *«). (3) 

The cancellation of acceptance effects has been verified using simulated events, including 
the full simulation of detector effects, trigger response and final event selection. Any 
non-cancelling acceptance bias on the measured lifetime is found to be smaller 1 fs. 

In order to extract the effective B® — > K + K~ lifetime, the yield of B meson candidates 
is determined in bins of decay time for both decay modes. Thirty bins between -1 ps and 
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35 ps are chosen such that each bin contains approximately the same number of B meson 
candidates. The ratio of the yields is then fitted as a function of decay time and the relative 
lifetime can be determined according to Eq. [3j With this approach it is not necessary 
to parametrise the decay time distribution of the background. In order to maximise 
the statistical precision, both steps of the analysis are combined in a simultaneous fit 
to the K + K~ and K + ix~ invariant mass spectra across all decay time bins. The signal 
distributions are described by Gaussian functions and the combinatorial background by 
first order polynomials. The parameters of the signal and background probability density 
functions (PDFs) are fixed to the results of time-integrated mass fits before the lifetime 
fit is performed. The B° — > K + tt~ yield (Ns^Kn) is allowed to float freely in each bin 
but the B® — >• K + K~ yield (N Bs ^ KK ) is constrained to follow 

N Bt -+KK(ti) = NB^tjRWe-^Kx ~ T ^), (4) 

where is the mean decay time in the i th bin. In total the simultaneous fit has 94 free 
parameters and tests using Toy Monte Carlo simulated data have found the fit to be 
unbiased to below 1 fs on the measured B® — > K + K~ lifetime. Each mass fit used in the 
simultaneous fit is unbinned and must be split into mass bins in order to evaluate the fit 
X 2 ■ Two mass bins are chosen, one signal dominated and one background dominated, in 
order to guarantee a minimum of 5-6 candidates in each bin. Using this appraoch the x 2 
per degree of freedom of the simultaneous fit is found to be 0.82. The right part of Fig. [I] 
shows the decay time distribution obtained from the fit and the fitted reciprocal lifetime 
difference is 

T KK ~ T kI = °- 013 ± °- 045 ( Stat ) PS^ 1 - 

Taking the B° — > K + 7t~ lifetime as equal to the mean B° lifetime (r^o = 1.519±0.007 ps) 
[8], this measurement can be expressed as 

t kk = 1.490 ± 0.100 (stat) ± 0.007 (input) ps 
where the second uncertainty originates from the uncertainty of the B° lifetime. 



4 Absolute Lifetime Measurement 

The absolute lifetime measurement method directly determines the effective B° s — > K + K~ 
lifetime using an acceptance correction calculated from the data. This method was first 



used at the NAll spectrometer at CERN SPS [12], further developed within CDF [13j[14 
and was subsequently studied and implemented in LHCb (l5fl6~] . The per event acceptance 
function is determined by evaluating whether the candidate would be selected for different 
values of the B meson candidate decay time. For example, for a B meson candidate, with 
given kinematic properties, the measured decay time of the B meson candidate is directly 
related to the point of closest approach of the final state particles to the associated primary 
vertex. Thus a selection requirement on this quantity directly translates into a discrete 
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decision about acceptance or rejection of a candidate as a function of its decay time. 
This is illustrated in Fig. [2] In the presence of several reconstructed primary interaction 
vertices, the meson may enter a decay-time region where one of the final state particles 
no longer fulfills the selection criteria with respect to another primary vertex. Hence the 
acceptance function is determined as a series of step changes. These turning points at 
which the candidates enter or leave the acceptance of a given primary vertex form the 
basis of extracting the per-event acceptance function in the data. The turning points are 
determined by moving the reconstructed primary vertex position of the event along the B 
meson momentum vector, and then reapplying the event selection criteria. The analysis 
presented in this paper only includes events with a single turning point. The drop of the 
acceptance to zero when the final state particles are so far downstream that one is outside 
the detector acceptance occurs only after many lifetimes and hence is safely neglected. 



(a) 
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Figure 2: Decay-time acceptance function for an event of a two-body hadronic decay. The 
light blue (shaded) regions show the bands for accepting the impact parameter of a track. 
The impact parameter of the negative track (IP2) is too small in (a) and lies within the 
accepted range in (b). The actual measured decay time lies in the accepted region. The 
acceptance intervals give conditional likelihoods used in the lifetime fit. 

The distributions of the turning points, combined with the decay-time distributions, 
are converted into an average acceptance function (see Fig. [3]). The average acceptance 
is not used in the lifetime fit, except in the determination of the background decay-time 
distribution. 

The effective B° s — > K + K~ lifetime is extracted by an unbinned maximum likelihood fit 
using an analytical probability density function (PDF) for the signal decay time and a non- 
parametric PDF for the combinatorial background, as described below. The measurement 
is factorised into two independent fits. 

A first fit is performed to the observed mass spectrum and used to determine the signal 
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and background probabilities of each event. Events with B° s candidates in the mass range 
5272 — 5800 MeV/ c 2 were used, hence reducing the contribution of partially reconstructed 
background and contamination of B° decays below the B° s mass peak. The signal distri- 
bution is modelled with a Gaussian, and the background with a linear distribution. The 
fitted mass value is compatible with the current world average (§]. 

The signal and background probabilities are used in the subsequent lifetime fit. The 
decay-time PDF of the signal is calculated analytically taking into account the per-event 
acceptance and the decay-time resolution. The decay-time PDF of the combinatorial 
background is estimated from data using a non-parametric method and is modelled by a 
sum of kernel functions which represent each candidate by a normalised Gaussian function 
centred at the measured decay time with a width proportional to an estimate of the density 
of candidates at this decay time 17 . The lifetime fit is performed in the decay-time range 
of 0.6 — 15 ps, hence only candidates within this range were accepted. The analysis was 
tested on the B° — > K + ir~ channel, for which a lifetime compatible with the world average 
value was obtained, and applied to the B° s — > K + K~ channel only once the full analysis 
procedure had been fixed. The result of the lifetime fit is 



r KK = 1.440 ± 0.096 (stat) ps 



and is illustrated in Fig. |3j 
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Figure 3: Left: Average decay-time acceptance function for signal events, where the error band is 
an estimate of the statistical uncertainty. The plot is scaled to 1 at large decay times, not taking 
into account the total signal efficiency. Right: Decay-time distribution of the B® — > K + K~ 
candidates and the fitted functions. The estimation of the background distribution is sensitive 
to fluctuations due to the limited statistics. Both plots are for the absolute lifetime measurement. 



5 Systematic Uncertainties 

The systematic uncertainties are listed in Table [T] and discussed below. The dominant 
contributions to the systematic uncertainty for the absolute lifetime measurement come 
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Table 1: Summary of systematic uncertainties on the B® — > K + K lifetime measurements. 



Source of uncertainty 


Uncertainty on 


Uncertainty on 




t kk (fs) 




Fit method 


3.2 




Acceptance correction 


6.3 


> 0.5 


Mass model 


1.9 


J 


B— > h + h'~ background 


1.9 


1.4 


Partially reconstructed background 


1.9 


1.1 


Combinatorial background 


1.5 


1.6 


Primary vertex association 


1.2 


0.5 


Detector length scale 


1.5 


0.7 


Production asymmetry 


1.4 


0.6 


Minimum accepted lifetime 


1.1 


N/A 


Total (added in quadrature) 


8.4 


2.7 


Effective lifetime interpretation 


2.8 


1.1 



from the treatment of the acceptance correction (6.3 fs) and the fitting procedure (3.2 fs). 
The systematic uncertainty from the acceptance correction is determined by applying the 
same analysis technique to a kinematically similar high statistics decay in the charm 
sector (D° — > K~n + |18|). This analysis yields a lifetime value in good agreement 
with the current world average and of better statistical accuracy. The uncertainty on the 
comparison between the measured value and the world average is rescaled by the B meson 
and charm meson lifetime ratio. The uncertainty due to the fitting procedure is evaluated 
using simplified simulations. A large number of pseudo-experiments are simulated and 
the pull of the fitted lifetimes compared to the input value to the fit is used to estimate 
the accuracy of the fit. These sources of uncertainty are not dominant in the relative 
method, and are estimated from simplified simulations which also include the systematic 
uncertainty of the mass model. Hence a common systematic uncertainty is assigned to 
these three sources. 

The effect of the contamination of other B — > h + h'~ modes to the signal modes is 
determined by a data-driven method. The misidentification probability of protons, pions 
and kaons is measured in data using the decays — > 7r + 7r~, D° — > K + n~ , (f> — > K + K~ 
and A — > pn~ , where the particle type is inferred from kinematic constraints alone [19] . As 
the particle identification likelihood separating protons, kaons and pions depends on kine- 
matic properties such as momentum, transverse momentum, and number of reconstructed 
primary interaction vertices, the sample is reweighted to reflect the different kinematic 
range of the final state particles in 5—^ h + h'~ decays. The effect on the measured lifetime 
is evaluated with simplified simulations. 

Decays of B° s and B° to three or more final state particles, which have been partially 
reconstructed, lie predominantly in the mass range below the B° s mass peak outside the 
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analysed region. Residual background from this source is estimated from data and eval- 
uated with a sample of fully simulated partially reconstructed decays. The effect on the 
fitted lifetime is then evaluated. 

In the absolute lifetime measurement, the combinatorial background of the decay time 
distribution is described by a non-parametric function, based on the observed events with 
masses above the B® meson signal region. The systematic uncertainty is evaluated by 
varying the region used for evaluating the combinatorial background. In the relative 
lifetime measurement, the combinatorial background in the hh' invariant mass spectrum 
is described by a first order polynomial. To estimate the systematic uncertainty, a sample 
of simulated events is obtained with a simplified simulation using an exponential function, 
and subsequently fitted with a first order polynomial. 

Events may contain several primary interactions and a reconstructed B meson can- 
didate may be associated to the wrong primary vertex. This effect is studied using the 
more abundant charm meson decays where the lifetime is measured separately for events 
with only one or any number of primary vertices and the observed variation is scaled to 
the B meson system. 

Particle decay times are measured from the distance between the primary vertex and 
secondary decay vertex in the silicon vertex detector. The systematic uncertainty from this 
source is determined by considering the potential error on the length scale of the detector 
from the mechanical survey, thermal expansion and the current alignment precision. 

The analysis assumes that B° s and B° s mesons are produced in equal quantities. The 
influence of a production asymmetry for B® mesons on the measured lifetime is found to 
be small. 



In the absolute lifetime method both a Gaussian and a Crystal Ball mass model 20 
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implemented and the effect on fully simulated data is evaluated to estimate the systematic 
uncertainty due to the modelling of the signal PDF. In the relative lifetime method this 
uncertainty is evaluated with simplified simulations and included in the fitting procedure 
uncertainty. 

In the absolute B® — > K + K~ lifetime measurement a cut is applied on the minimal 
reconstructed decay time. As the background decay time estimation will smear this step 
in the distribution, a systematic uncertainty is quoted from varying this cut. 

There is an additional uncertainty introduced if the result is interpreted using Eq. [2j as 
this expression does not take into account detector resolution and decay time acceptance. 
This effect was studied using simplified simulations modelling the acceptance observed 
in the data and conservative values of Ar s = 0.1 ps and Aat s = -0.6. The observed 
bias with respect to the prediction of Eq. [2] is 3 fs. This effect is labelled "Effective 
lifetime interpretation" in Table [T] and is not a source of systematic uncertainty on the 
measurement but is relevant to the interpretation of the measured lifetime. 
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6 Results and Conclusions 



The effective B® — > K + K~ lifetime has been measured in pp interactions using a data 
sample corresponding to an integrated luminosity of 37 pb -1 recorded by the LHCb ex- 
periment in 2010. Two complementary approaches have been followed to compensate 
for acceptance effects introduced by the trigger and final event selection used to enrich 
the sample of B® mesons. The absolute measurement extracts the per event acceptance 
function directly from the data and finds: 

Tkk = 1-440 ± 0.096 (stat) ± 0.008 (syst) ± 0.003 (model) ps 

where the third source of uncertainty labelled "model" is related to the interpretation of 
the effective lifetime. 

The relative method exploits the fact that the kinematic properties of the various 
B— > h + h'~ modes are almost identical and extracts the B° s — > K + K~ lifetime relative to 
the B°—> K + 7t~ lifetime as: 

t kk - t kI = °- 013 ± °-°45 (stat) ± 0.003 (syst) ± 0.001 (model) ps" 1 . 

Taking the B°— >■ K + tt~ lifetime as equal to the mean B° lifetime (r^o = 1.519±0.007 ps) 
[8], this measurement can be expressed as: 

tkk = 1-490 ± 0.100 (stat) ± 0.006 (syst) ± 0.002 (model) ± 0.007 (input) ps. 

where the last uncertainty originates from the uncertainty of the B° lifetime. Both mea- 
surements are found to be compatible with each other, taking the overlap in the data 
analysed into account. 

Due to the large overlap of the data analysed by the two methods and the high 
correlation of the systematic uncertainties, there is no significant gain from a combination 
of the two numbers. Instead, the result obtained using the absolute lifetime method is 
taken as the final result. The measured effective B° s — > K + K~ lifetime is in agreement 
with the Standard Model prediction of tkk = 1-390 ± 0.032 ps [3]. 
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